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Abstract 

The effective action for wee partons in large nuclei includes a sum 
over static color sources distributed in a wide range of representa- 
tions of the SU(N C ) color group. The problem can be formulated as 
/\ • a random walk of partons in the N c — 1 dimensional space spanned 

by the Casimirs of SU(N C ). For a large number of sources, k 3> 1, 
we show explicitly that the most likely representation is a classical 
representation of order O(vfc). The quantum sum over representa- 
tions is well approximated by a path integral over classical sources 
with an exponential weight whose argument is the quadratic Casimir 
operator of the group. The contributions of the higher N c — 2 Casimir 
operators are suppressed by powers of k. Other applications of the 
techniques developed here are discussed briefly. 



1 Introduction 

In QCD, there are a number of physical situations where higher dimensional 
representations of the gauge group are relevant. In particular, in perturba- 
tive QCD, because of the ubiquitous nature of bremsstrahlung processes, one 
often encounters situations where partons at one energy scale interact simul- 
taneously with several partons at a different energy scale [1]. The former 
may therefore couple to color charges in a wide range of representations of 
the color group. If one can argue that higher dimensional representations are 
the most likely representations, it may be possible, by the correspondence 
principle, to treat these as classical representations. Classical color charges 
are considerably simpler to treat than their quantum counterparts. 

In this paper, we will address the following formal questions which are 
relevant to addressing this issue quantitatively. Given k non-interacting 
quarks in the fundamental SU(N C ) representation, what is the distribution 
of degenerate irreducible representations that one generates? What is the 
most likely representation? Is this representation a classical representation? 
It is well known that the N c — *■ oo limit of QCD is a classical theory even for 
k — 1 [2]. What happens for finite N c but k ^> 1 x ? We will demonstrate in 
the following that, for an SU(N C ) gauge theory, the problem can be formu- 
lated as a random walk problem in the space spanned by the N c — 1 Casimir 
operators of the group. 



lr rhe former corresponds to the large N limit where the invariance group of the theory 
grows with N. Several such large N cases have been discussed in Ref. [2]. We are instead 
interested in a situation where the underlying symmetry group is unchanged but where 
quantum operators in higher dimensional representations appear. Examples of these are 
the large spin limits of quantum spin models [3] . 



Our interest in the questions posed here arise in the context of a par- 
ticular model of small x physics in QCD, the McLerran-Venugopalan (MV) 
model for small x parton distributions in large nuclei [4]. The MV model 
is a simple model to understand the physics of saturation [5] at small x in 
QCD. To highlight the physical relevance of the formal mathematics, we will 
therefore address solutions to these general questions in the specific context 
of the MV model. However, the solution to the mathematical problem posed 
may be of more general interest and will hopefully prove useful in developing 
novel treatments for a variety of physical situations in QCD. These include 
problems in jet physics [6], finite temperature transport theory [7] and even 
percolation and string based models of multi-particle production [8, 9]. 

This paper is organized as follows. In the following section, we will 
briefly review the McLerran-Venugopalan model. We will focus in particular 
on those aspects of the model related to the random walk of partons in the 
space spanned by the Casimirs of their color group and of their subsequent 
treatment as classical color charge distributions. In section 3, we consider 
the case of the MV-model in an SU(2) gauge theory. The problem here 
is especially simple since it can be mapped into a one dimensional random 
walk problem of k spins of spin-1/2. The peak of the distribution is at 
a representation whose dimension is of order 0(\/k). The distribution of 
representations about this peak value is well approximated as an exponential 
in the quadratic Casimir operator weighted by k, as indeed assumed in the 
MV-model. 

The "real world" case of an SU(3) gauge theory is discussed in section 
4. In this case, one can think of the random walk as occurring in the space 



spanned by the third component of the isotopic spin J 3 and the hypercharge 
Y. Determining the distribution of higher dimensional representations gen- 
erated by k quarks when the random walk is in more than one dimension 
(SU(N C ) for iV c > 3) is non-trivial and is the primary subject of this paper. 
We will show explicitly, that for a random walk of k quarks in the fundamen- 
tal SU(3) representation, the distribution of degenerate representations is, 
as in the £77(2) case, also peaked at a representation of order 0{yk). The 
distribution of representations about the most likely representation, how- 
ever, cannot simply be expressed in terms of an exponential of the quadratic 
Casimir operator alone; the argument of the exponential has an additional 
term proportional to the cubic Casimir operator weighted by k 2 . The cubic 
term thus acts as a perturbation of the quadratic term. The relative magni- 
tude of the former to the latter is of order 0(1/ yk) with a coefficient that is 
computed exactly. Thus for very large nuclei A — ► oo, the contribution from 
this term can be neglected 2 . We also consider the representations generated 
by random walks of gluons and of quark- ant i- quark pairs. Interestingly, in 
both of these cases, the distribution of representations is given by an expo- 
nential in the quadratic Casimir alone. This is true because the quadratic 
Casimir, unlike the cubic Casimir, is symmetric in the upper and lower £77(3) 
tensor indices (denoted by m and n in the paper). Gluons and quark-anti- 
quark representations, unlike the quark representations, favor representations 
that are symmetric in m and n. 

The generalization of these results to SU(N C ) for any general N c is 

straightforward and is briefly discussed in section 5. The essential features of 

2 The statements here are valid for values of x where small x quantum evolution is not 
too important. Quantum evolution effects will be discussed briefly later in this section 
and in the final section. 



the SU(3) result are preserved: the distribution of representations about the 
most likely representation is determined primarily by the quadratic Casimir. 
The other N c — 2 Casimirs add small (and in principle quantifiable) pertur- 
bative corrections-all of which vanish for k — ■> oo. 

In section 6, we will summarize our results and discuss their ramifica- 
tions. A brief discussion of classical limits of quantum systems in general, 
and quantum spin systems in particular, is presented in appendix A. The 
other appendices contain technical details of results presented in the body of 
the paper. 

2 Classical color charges in the McLerran- 
Venugopalan model 

The McLerran-Venugopalan model is a classical effective field theory for wee 
parton distributions in large nuclei [4]. The model is formulated in the light 
cone gauge A + = 0, and in the infinite momentum frame, where the mo- 
mentum of the nucleus, P + — ► oo. In this case, the physics of time dilation 
ensures that the time scales for partons carrying a higher fraction x of the 
nuclear momentum (the "valence" partons with x ~ 1) to interact with one 
another is much larger than the characteristic time scale for their interac- 
tions with softer "wee" partons (with i< 1). The softer partons, in turn, 
couple to a large number of valence partonic sources, which appear static 
on the wee parton time scales. In particular, wee partons with coherence 
lengths (/ co h. ~ l/^m^x) much greater than the Lorentz contracted nuclear 
widths (of order 2R/ n f), or equivalently with x <C A^ 1 ^ 3 , couple coherently 
to ~ A 1 / 3 valence partons along the longitudinal length of the nucleus in the 



infinite momentum frame. Here m^ is the nucleon mass and 7 ~ P + /m^ is 
the Lorentz factor in the infinite momentum frame. 

The separation of scales in x between wee and valence partons is experi- 
mentally well established in parton distribution function measurements [10]. 
The insensitivity of the valence distributions to the sea is also experimentally 
well established-the well known phenomenon of "limiting fragmentation" is 
a direct consequence [11]. 

Further, on the time scales relevant to wee parton dynamics, it is reason- 
able to assume that the valence parton sources are random light cone sources. 
This is plausible because firstly, most of the hard partons are confined in dif- 
ferent nucleons (and hence do not interact with one another); secondly, due 
to time dilation, even hard partons in the same nucleon are non-interacting 
(and therefore independent sources) over the short time scales relevant to the 
wee partons. 

How many of these random sources the wee partons actually couple to 
depends on the typical transverse momentum of the wee parton 3 . A wee 
parton with momentum p± resolves an area in the transverse plane (A:r_i_) 2 ~ 
l/p±- The number of valence partons it interacts simultaneously with is then 

k = k (Ax±) 2 = ^g cc (Ax ± ) 2 , (1) 

which indeed is proportional to A 1 / 3 since N va \ ence = 3 • A in QCD. This 
counting of color charges is only valid as long as p± > Aqcd since, on the scale 
of the nucleon size , p± ~ Aqcd ~ 200 GeV, confinement ensures that the wee 
partons see no net color charge. The momentum scale for color neutrality may 



3 The wee parton is soft only in its longitudinal momentum-its transverse momentum 
may be large. 
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be even larger than the confinement scale due to screening effects analogous 
to Debye screening [12, 13]. These effects, while very interesting, are beyond 
the scope of this paper. 

It was argued in the MV model that when k ~ A 1 ! 3 ^> 1, the most likely 
color representation that the wee partons couple to is a higher dimensional 
representation. We will demonstrate later that this representation is one of 
order 0(\/k). Thus, as discussed further in appendix A, for large enough 
k, the color charge distribution can be treated as a classical color charge 
distribution. It is further argued in the MV-model, that the distribution of 
classical representations is Gaussian, with a variance proportional to k (or, 
equivalently, to A 1 ' 3 ). 

With these stated assumptions, the classical effective Lagrangian for the 
MV-model, formulated in the infinite momentum frame (P + — »■ oo) and light 
cone gauge (A + = 0) has the form 



C = <fx 



-F a F a -J -A 
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2/4' 

(2) 



The first term is the usual QCD Field Strength tensor squared (Lorentz 
indices are suppressed for simplicity here), which describes the dynamics of 
the wee partons. The second term denotes the coupling of these wee partons 
to the hard valence parton sources. The valence parton current has the form 

where p is the classical color charge per unit transverse area of the valence 

parton sources 4 . The likelihood of configurations of differing p's is given by 

4 As written, this term in the action is not gauge invariant. Gauge invariant expressions 
are given in Refs. [14, 16] . The lowest order in g term in the expansion of these expressions 
gives the form shown here. 



the final term in Eq. (2) - the weight n\ of the Gaussian is the average color 
charge squared per unit area per color degree of freedom. The correlator of 
color charge densities is then 5 

(p a (x±)p\y±)) = p 2 A 5 ab 5^(x ± - y ± ) . (4) 

The average color charge squared per unit area per color degree of free- 
dom, n\ is simply determined in the MV-model. The color charge squared in 
a tube of transverse area (Axj_) 2 is the color charge squared per unit quark 
(g 2 Cp) times the number of quarks in the tube: (A:r_|_) 2 • AN c /ttR 2 . When 
this is normalized as in Eq. (4), per unit transverse area, per color degree of 
freedom, one obtains 

2 9 2 A , , 

^=2^' (5) 

which is of order A 1 / 3 . For a very large nucleus, fx\ ^> Aq CD , and since 
it is the only scale in the problem, the coupling constant a s = a s (fi\) <C 
1. Since the coupling is weak, one can compute parton distributions for a 
large nucleus. The classical field equations can be solved and it was shown 
explicitly that the number distribution is of order 1/as for p\ < g 2 fi\ [18]. 
The Gaussian functional weights and the classical field equations for a 
large nucleus are also recovered in a particular model of large nuclei where 
the nucleons are modeled as color singlet quark anti-quark pairs [19]. In 
Refs. [18] and [19], it was recognized that it was essential to smear the 
(^-function sources in Eq. (3) in the x~ direction to obtain regular classical 

solutions. Quantum corrections to the MV-model are large [20] but they can 

5 There are several conventions for the color charge densities in the literature. For a 
discussion, sec Rcf. [171. 



be absorbed, via a Wilsonian renormalization group procedure, at each step 
in x into new sources and fields, while preserving the essential structure of 
the classical effective Lagrangian [14, 21]. Thus this simple MV-model can be 
reformulated as a more sophisticated effective field theory, the "Color Glass 
Condensate" (CGC), which is applicable to both hadrons and nuclei. In the 
CGC, the weight functional can be represented more generally as 

exp(-W\p]), (6) 

where W[p] is a gauge invariant functional that satisfies the non-linear renor- 
malization group (RG) equation in x. For a recent review of the MV model 
and the CGC, see Ref. [22]. 

In the rest of the paper, we will restrict ourselves to the MV model and 
seek to establish more rigorously the following, namely, 

• the most likely representation the wee partons couple to is a higher 
dimensional representation which can be represented in terms of a clas- 
sical color charge density p, 

• the sum over color charges in the path integral can be expressed as 
a path integral over classical color charges, with a weight that is well 
approximated by a Gaussian for a large number k of color charges. This 
is shown for quark sources as well as for gluon and quark-anti-quark 
sources for N c > 3. 

• Finally, we will discuss, for N c > 3, possible sub- leading corrections to 
the Gaussian term and their relative dependence on k. 



3 Random walks and classical color charge 
representations in two color QCD at small 

x 

We now consider the case of the MV-model in SU(2) QCD. As discussed 
previously, (see Eq. (1)) the wee partons couple to a large number of uncor- 
rected quarks-in this case, quarks in the fundamental SU(2) representation. 
We wish to discuss here the distribution of representations generated by 
adding A;-quarks. 

3.1 Random walk of many spin- 1/2 quarks 

The problem of adding k random color charges in N c = 2 QCD is exactly 
equivalent to the problem of adding k spins of spin-l/2-they both correspond 
to an internal symmetry group whose generators are elements of the SU(2) 
algebra. Therefore, let's start with a spin singlet state. Multiplying with 
a spin 1/2 state results in a 1/2 state, or 

2x1 = 2 (7) 

where 1 is the singlet state and 2 denotes the spin 1/2 state. In general we 
denote a spin I state by the degrees of freedom associated with the state, 
namely, 8 — 21 + 1. 

Continuing, multiplying another spin 1/2 states results in 

2x2 = 1 + 3 (8) 

If we have 3 spin 1/2's, we get 

2x2x2 = 2x (1 + 3) = 2 + 2 + 4, (9) 

10 



and so on. 

Let v^ denote the multiplicity of the representation s when k funda- 
mental representations are multiplied. In general, if one adds one more spin 
1/2 particle to an s state, the result is a mixture of the states s — 1 and 
s+1: 

2 x s= (a- l) + + l) (10) 

Therefore, the multiplicity of s state in the k-th iteration is given by 

vP=v^+v^ (11) 

We wish to find the solution of this recursion relation with the initial 
condition 

4 0) = 1 otherwise ^ 0) = (12) 

We begin by noting that the binomial coefficients satisfy a similar (but not 
identical) recursion relation, 

(D-fciM*: 1 )- 

Therefore, we can try and find a solution to Eq. (11) in terms of binomial 
coefficients. 
Define 

G «"=((t+*.)/2)- (M) 

which is non-zero only if (k + s) is even. Consider now the combination 

g(k : s) = Gk-.s-x -G k:s+1 

= { (k + s -l)/2 J " { (k + s + l)/2 J ■ (15) 

11 



Using the binomial coefficient relation (13), one can easily show that this 
combination satisfies the recursion relation 

g(k:s)=g(k-l:s-l)+g{k-l:s + l) (16) 

Furthermore, for k = 

and all the other g(0 : s)'s are zero. 

In this way, we have found a solution to the recursion relation (11) 
satisfying the right boundary condition, Eq. (12), namely, 

v^=g(k:s) = G k:s ^-G k .. s+1 (18) 

Since we know the G's explicitly in terms of binomial coefficients (Eq. (14)), 

we can determine v^ explicitly. Since we are interested in the limit k ^> 

s ^> 1, we can use the well known Stirling formula to approximate 

k\ 
Gk - = ((k + s)/2)\((k-s)/2)\> (19) 

as 

lnG fc:s « A; In 2 - s 2 /2k - (1/2) In A; - (1/2) ln(27r) , (20) 

or equivalently, 

ofe-l/2 

G k .. s ^^=e~ s2 l 2k . (21) 

vkn 



Thus the solution of the recursion relation for k ^> s ^> 1 is 

Vg = Gk-s-i — Gk-s+i 
~ —2d s Gk: S 

nk+l/2 

se~ s2 l 2k (22) 



kykn 
12 



This formula represents the distribution of representations with spin I = 
(s — l)/2 when there are k random spin 1/2 particles in the system. 

To find the probability of the spin state I, we need to take into account 
the degeneracy of the spin state s = 2/ + 1. Since we have k random spin 
1/2 particles, summing the degeneracy times the multiplicity (st)f) over all 
possible representations must be equal to the 2 k degrees of freedom. Even 
with our large k approximation, we find that this is indeed the case: 

dssvf ) = 2 k . (23) 

o 

Defining the "classical" spin vector 6 1 which is related to s by the 
relation 

s 2 = 41 2 , (24) 

or equivalently l 2 = (/ + 1/2) 2 , we can exploit the lack of angular dependence 
in the argument to re-write Eq. (23) as 

/ 9 \ *^/^ r 
l=f— J / d 3 / exp (-2 1 2 /A;) . (25) 

Thus one can define 

/ o n, 3/2 

to be the probability density for the system of k spin- 1/2 particles to have 
the total spin 1. The above formula clearly resembles the classical Maxwell- 
Boltzmann distribution of particles in a heat bath. 



6 We will see a posteriori that 1 is of order k. For large k, we argue in Appendix A 
that such representations are classical representations. 
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Our result can be easily expressed in terms of the SU(2) Casimir, which 

is 

D 2 = 1(1 + 1) « l 2 . (27) 

The most probable value of the Casimir is then given by 

D 2 = \ (28) 

orl = Vk/2 + 0(l). 

One can repeat the same analysis for the case of many spin-1 particles. 
In this case, one has a somewhat more complex analysis (involving trinomial 
coefficients) but the final result for the multiplicity distribution is identical 
to the spin-1/2 case-up to overall constants. 

What does this result for v ^ imply for the derivation of the effective 
action in the MV model? In order to understand this, we have to step back a 
little and discuss, with some greater detail, the derivation of the MV-effective 
action [4]. 

3.2 From random walks to path integrals 

The path integral describing the ground state \0) of a large nucleus can be 
written as 



Z= (0\e lx " f qcd|0)= lim ^(7V,Q|e" P *™\N,Q) , (29) 



The sums over N and Q here respectively represent the sum over all possible 
states in the path integral and a sum over the color quantum numbers of 
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these states 7 . In the lattice representation of the path integral, a state 
\N,Q) would correspond to a 4-dimensional box on the lattice containing a 
net color charge Q. In the MV-model, one is interested in a coarse grained 
theory, where the size of the box is chosen (see the first paper in Ref. [4]) 
such that it contains a large number of quark color charges, namely, k >> 1. 
From Eq. (1), this implies that the coarse grained state \N, Q) contains only 
modes with Aq CD <C p\ <C fj,\ (see Eq. (5)) and the corresponding color 
charge Q is the charge constructed from all the k(A x± y 2 quarks contained in 
a box with transverse size l//x^ <C (Ax±) 2 <C 1/Aq CD . 

In general, one has to perform the quantum mechanical sum over the 
color charges to construct \N,Q), which is a difficult problem indeed. It 
is at this point one sees the relevance of the problem of the distribution of 
representations that we have just solved (for the N c = 2 case). In the MV 
model, the color charges are random, and k charges can be distributed in 
a wide range of distributions. Since the most likely representation is s = 
yk/2 ^> 1 for k >> 1, the color charge corresponding to this representation, 
as argued in appendix A, is a classical representation. 

Thus in the SU(2) case, the sum over all spin states in the path integral 
can be replaced by the integral 8 

5>f fc) £ \l,m)(l,m\^jdHe- 2l2 / k , (30) 

l m=—l 

where we have made use of Eqs. (22) and (26). For a transverse area of 

size (Ax_i_) 2 containing k spins, then one can introduce classical color charge 

7 The distinction between the quantum numbers represented by N and the color degrees 
of freedom represented by Q is of course artificial. We make this distinction because it 
will prove useful in constructing the effective theory. 

8 This assumes that the s = 21 + 1 states in a representation have the same P + . 

15 



(spin) density p a which is defined by the relation, 
1 r (Ax_l) 2 

l a = ~ / d 2 X ± p a (x ± ) W Pa(x±) • (31) 

9 J (^±) 2 9 

We can now re-express the Gaussian weight in Eq. (30) in terms of the 
classical color charge density: 

l 2 l 2 

= N c — — p a p a 

9 & 

= -^J^fPaPa- (32) 

9 A 
Here we have used Eq. (1) to express k in terms of the number of valence 
quarks in an area (Axj_) 2 . We also used the fact that the number of valence 
quarks per baryon is equal to N c . 

The sum in Eq. (30) can therefore be expressed as 

/ J] d P a W[p] ee J n dp* exp (- f A d*x ± PiX± ^ X±) > j (33) 

where 



9 2 A 
2irR 2 



V-a = ±^ (34) 



is independent of iV c , and is the color charge squared per unit area defined 
in the MV-model (Eq. (5)). Integrating over all the boxes of area (A:r_|_) 2 in 
the transverse plane, one obtains the path integral for classical color charge 
distributions in the MV-model [4]. 
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4 Random walks and classical color charge 
representations for N c = 3 QCD at small x 

We shall now extend our discussion in the previous section to the case of 
three color QCD. The problem is more non-trivial here since one now has a 2- 
dimensional random walk in the space generated by the two SU(3) Casimirs. 
Thus it is no longer obvious that the quadratic Casimir alone will contribute 
to the "Boltzmann" weight in the path integral over color charges. We will 
show in this section that the approximation that only the quadratic Casimir 
contributes is an excellent approximation when the number of quarks is large. 
We will first consider the case where the sources of color charge are £77(3) 
quarks only before considering the representations generated by random glu- 
ons and quark-anti-quark states respectively. 

4.1 Many 577(3) quarks 

In the SU(3) case, representations are labeled by two integers m and n 
and written as (m, n). Here m and n are the number of the upper and 
the tensor indices, respectively. A fundamental quark 3 state in SU(3) is 
hence (1,0) while an anti-quark in the 3 state is (0,1). As is well known, 
SU(3) representations can be conveniently represented graphically as Young 
tableaux, where m denotes the number of boxes in the uppermost row minus 
those in the middle row, while n is the number of boxes in the middle row 
minus that of the bottom row. 

Recall that we are interested in studying the distribution of representa- 
tions when one adds k quarks in the fundamental representation. Consider 
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Figure 1: Young tableaux representing the result of multiplying SU(3) quarks 
in a representation labeled by the tensor indices m & n with a quark in the 
3 representation. 



what happens when one adds a quark to an (m, n) representation. This is 
represented by the Young tableaux in Fig. 1 which allows us to deduce that 

(1,0) X (m, n) = (m + 1, n) + (ra — 1, n + 1) + {m, n — 1) (35) 

We shall represent the multiplicity of each state by a matrix element 
N m} n with the indices < m, n representing the same (m, n) of the represen- 
tations. (Note that these N myn are the SU(3) analogs of the v s in the SU(2) 
case.) The multiplicity of the (m,n) state in the (k + l)-th iteration is then 
given by the recursion relation 



ju(fc+l) = lYf( k ) , lYf( k ) , N W 

m,n m— l,n ' 7n+l,n— 1 ' ly m,n+l 



(36) 



m 



rW 



for m, n > with the understanding that N^Li = N_{ n = 0. The initial 
condition is 



N t 



(0) 



0,0 



1 otherwise iV"W n = 



(37) 



As in the 577(2) case, we now wish to determine N^ n by solving this recursion 
relation. 

We begin by noting that Eq. (36) represents a trifucation process since 
we are multiplying 3's. The basic unit should therefore be the tri-nomial 
coefficients 9 : 

k\ 
III 111 h- 

Following the pattern set by Eq. (14) and Eq. (35), we define 

(to, n) = h(l, 0) + Z 2 (-l, 1) + / 3 (0, -1) , (39) 

Re-defining Ck;h,i 2 ,h * n terms of to, n as Gk- m , n , we have, 

( k+2m+n \ | / k—m+n \ | / h—m—ln \\ 
\ 3 )■{ 3 )■{ 3 )■ 

We are now in a position to compute the multiplicity of representations. One 
finds that 

m,n ~ ^fe: m,n r t-Jfc: m+3,n "T >-JA;: m,n+3 

— Gk: m+2,n-l ~ Gk: rn-l,n+2 ~ Gk: m+2,n+2 (41) 

satisfies the recursion relation Eq. (36) and the initial condition Eq. (37). We 

have verified this analytically and have checked numerically as well that the 

recursion relations (36) and our solution (41) indeed yield the same values 

upon iteration. 

As performed previously in the SU(2) case, we can use the Stirling 

formula to write the above expression for G (after some algebra) compactly 
9 Trinomial coefficients also arise in the adjoint spin-1 representation of SU(2). 
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as 

/3 fc+3 / 2 \ (m 2 + mn + n 2 ) (m — n)(2m + n){m + 2n) 
In G fe; m>B » In \^- j + — (42) 

for k >> m,n ^> 1. 

Recall that the S't/(3) Casimir for a given m, n state is given by [23] 

D (m,n) = W m 2 + mn + ^ + ^ + ^ _ ^3^ 

Note further that this Casimir is symmetric in m and n. This is unlike the 
Cubic Casimir which is ant i- symmetric under exchange of m and n and can 
be written as [24] 

D (ra,n) = _ ( m + 2n + 3) ( n + 2m + 3) (m - n) (44) 

18 
A little algebra should suffice to convince ourselves that G can be simply 
written in terms of the SU(3) Casimirs as 10 

G k:m , n « ^ exp (-3D" 2 l ' n /k) (l + 3D™' n /k 2 ) . (45) 

Clearly, near the peak, one has m = 0(yk) and n = 0(yk). Hence the 
cubic Casimir introduces a correction of size 0(l/\/k) for large k. 

In the large k limit, one finds that the multiplicity can be approximated 

as 

m n ^k: m,n t *J"fc: m+3,n > >-Jfc: m,n+3 

~~ Gk: m+2,n-l ~ Gk- m-l,n+2 ~ Gk: m+2,n+2 

27 mn (m + n) 3 2 +k qnm ,„ /n ,„ R x 
k s 2~I^ 6XP ( 2 ,k) ' ( } 



10 Here we have kept only the leading terms in D2 and D3 and expanded the exponential 
in D3 to lowest order. Our approximation can be checked to be self-consistent. 
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where we have dropped the D 3 contribution. 

Recall that the dimension of an SU(3) representation is 

(m + l)(n + l)(m + n + 2) mn(m + n) , . 

Therefore, just as in the SU(2) case, the probability to find the system of 

k quarks in an (m,n) representation should be proportional to d mn Nffi n . 

Furthermore, one should recover the total number of degrees of freedom 3 h 

by summing d mn Nffl n over all possible representations. 

With the expression obtained so far, one can show that integrating NJ£\ 

weighted by d mn over < m, n < oo yields, 

97 x 3§ +k r°° 

—^ / dmdnm 2 n 2 (m + nf exp (-3D™' n /k) = 3 k+1 (48) 

4 K 4 IX Jo 

The discrepancy with the integral above arises from the fact that the 
actual N^' n is rather asymmetric in m and n n . Since we are multiplying k 
number of 3 = (1,0) representations, the m side is more heavily populated 
than the n side. On the other hand, the approximation given in Eq. (46) is 
symmetric in m and n which is valid only in the vicinity of the peak. 

Nevertheless, our result enables us to normalize the distribution of rep- 
resentations as 

1 = — — — / dmdnm n (m + ri) exp (— ?>D™ ,n /k) (49) 

4 At 7i Jo 

We shall show below that we can re- write Eq. (49), in analogy to 

Eq. (25), as 

1=(^) 4 jd 8 Qe- N ^ 2 / k (50) 



11 Our method yields the correct number of degrees of freedom for gluons and quark- 
anti-quark pairs which, unlike quarks, are symmetric in m and n (see sections 4.2 and 
4.3). 
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where Q = (Q±, Q 2 , • ■ • , Qs) is a classical color charge vector defined by 



|Q| = \JQ a Q a = v/ D™ ,n and Qi,- ■ ■ ,Qs are its eight components. Once we 
are able to do this, the rest of the argument, expressing the sum over color 
charges in large nuclei as a path integral over Gaussian color charges, will 
follow in smooth analogy to the derivation in section 3.2. 

For any SU(3) representation (not necessarily a classical representa- 
tion), one can write the measure d 8 Q as [25, 26, 27] 

d 8 Q = dfadfadfaditid^ditsdm dn mn(m + n)— — , (51) 

V 48 / 

where fa and i\i (the so called "Darboux" variables) are canonically conjugate 

variables 12 . If the integrand depends only on m and n, we can carry out the 

integral over the angles dfa and their canonically conjugate momenta diii. 

One obtains (see appendix B for details of the derivation) [26] 

r 3 (2tt) 3 

/ II °^ ^ 7r * = mn{m + n) . (52) 

With this result, we can write Eq. (51) as 

d 8 Q = i-^= f dm dn (m 2 n 2 {m + nf) . (53) 

Substituting Eq. (53) into Eq. (49), we obtain Eq. (50). 

In summary, we have shown in this section that for a sum over classical 
SU(3) representations, the measure can be expressed as an eight dimensional 
integral over the components of the classical color charge Q (with a magni- 
tude of order yk ^> 1). 



12 It can be checked that Q\ ■ ■ ■ Qs satisfy {Q a , Qb} PB — fabcQc where the Poisson 
Bracket is defined as {Q a , Q b } PB = £? (fgf §& - f§f f 
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4.1.1 From random walks to path integrals in SU(3) QCD 

The reader may note that Eq. (50) is exactly analogous to Eq. (25) obtained 
for SU(2) quark sources. Then replacing 1 — > Q in Eqs.(30), (31) and (32), 
we recover 

/ J] dp a W[p] ee J II d P a exp (- J a d 2 x ± P(X± ^ X±) ) (54) 



where now a = 


= !,-• 


• , 8 and 


Va = 


9 2 A 
2itRY 



(55) 

is the color charge squared per unit area and is independent of N c as stated 
in the previous section. Thus, even in the SU(3) case, and as conjectured 
in the MV model, the leading contribution to the path integral measure 
over classical color charge densities has a Gaussian weight proportional to 
the quadratic Casimir. As shown in Eq. (45), the contribution of the cubic 
Casimir to the weight is 0(l/yk) which vanishes when k — ► oo. 

4.2 Many 577(3) gluons 

In the case of uncorrelated gluons, we are multiplying the 8's. Since the 
Young tableaux for one 8 has three boxes, we add three boxes at a time as 
we build up higher dimensional representations. As a consequence, there are 
sites in the (to, n) plane we will never visit in the random walk. One can see 
this as follows. In general, a Young tableaux diagram for a state (to, n) will 
have to + 2n + 3p boxes, where p is an integer with p > 0. Therefore, since 
we are multiplying boxes in multiples of 3, unless to + 2n = 3/ (where again, 
I is an integer > 0), the site (to, n) can never be reached from the origin. In 
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general, the allowed states in the (to, n) plane have the form (3s + u,3t + u), 
where s, t are integers > 0, and < u < 2. Thus for example, the state (1, 1) 
is allowed but not (1,2) or (1,0). 
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Figure 2: Young tableaux denoting the result 
of the multiplication of an arbitrary SU(3) 
representation labeled by (m,n) with an 8 
representation 



Fig. 2 shows Young tableaux illustrating the representations generated 
by multiplying an arbitrary 577(3) representation (to, n) with a (1, 1) (or 8) 
representation. It is clear from the multiplication of Young tableaux that 
one obtains 

(1, 1) X (m,n) 
= (m + l,n+ 1) + (ra + 2,n - 1) 
+ (m + 1, n - 2) + (m - 1, n + 2) 
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+ (m, n) + (m, n) + (m — 2, n + 1) + (wi — 1, n — 1) 

(56) 

Just as in the case of many quarks we dealt with previously, we will denote 
the multiplicity of a given state (m, n) by a matrix N myn with m, n > 0. Our 
initial state matrix is 

N$ = 1 otherwise N$ n = (57) 

From Eq. (56), one can deduce immediately that the multiplicity of states is 
given by recursion relation, 

Artfc+i) = 2 /V (fc) 4- /V (fc) 4- /V (fc) 

m,n m,n ' *m+2,n— 1 ' 1 *m+l,n+l 

+ Jv m+l,n-2 + iv m-l,n+2 + iv m-l,n-l + iv m-2,n+l l oc V 

Note that only states with m and n > are allowed, so in the above, the 
multiplicity of states whose labels are negative is zero. 

There are a few special cases which do not satisfy the above recursion 
relation. These are the states that lie on the edges of the first quadrant, 
namely, (m > 0,0) and (0,n > 0) and at the origin (0,0). These satisfy the 
following relations: 

< +1) = <i } , (59) 

This is because when repeatedly multiplying the 8's, a singlet state can only 
come from 8x8. For m > and n > 0, 

<r } = N% + ivi fc | M + N^_ 1>2 + N^_ 2tl (60) 

< +1) = <2 + iVffi-i + < } + i + <t 2 (61) 
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From Eq. (58), one can guess that the basic unit in solving the recursion 
relation should be 

Ck - hM '- M " yy-y(fc-ELy! 2 " Etl * ' (62) 

which appears in the expansion of 8 fe = (1 + 1 + 1 + 1 + 1 + 1 + 2) k . 
Following the pattern established thus far in Eq. (58), we set 

(m, n) = h(l, 1) + Z 2 (-l, -1) + l 3 (2, -1) + / 4 (-2, 1) + l 5 (l, -2) + Z 6 (-l, 2) 

(63) 

One can then define the following function (again in analogy to section 4.1), 

k 

Gk: m,n = /2 ^k: h,h,-,l 6 ^m-[(h-l 2 )+2(l 3 -l 4 )+(l 5 -l 6 )] ^n-[(h-l 2 )-(l :i -h)-2(h-l fi )] 
h,-,k=0 

(64) 

The solution to Eqs.(58), (59), and (61) is again given by the same combi- 
nation of Gfc:m,ra's as in the quark case 

m n ^k: m,n ~T~ Lrfc : m+3,n "T Lrfc : m,n+3 

~~ Gk: m+2,n-l ~ Gk: m-l,n+2 ~ Gk: m+2,n+2 ■ (65) 

The validity of this form of solution has been checked extensively by numer- 
ical means. Namely, we generated NJ£\ numerically by using the recursion 
relations and compared it with the values obtained by the solution (65). We 
are therefore confident that Eq. (65) is the solution. 

To express Gk-. m ,n i n terms of the SU(3) Casimirs, we can use the fol- 
lowing integral representation of the Kronecker delta: 

5„,= r^^-\ (66) 

26 



which enables us to carry out the sum in Eq. (64). This yields 

_ r* d$ r* d<p 

k:m ' n i-.27ri-.27r 

x e i<t>m+ilpn [2 + 2 cos(0 + <p) + 2 cos(20 -<p) + 2 cos(2y9 - <j>) J (67) 
Using the small angle approximation (valid for large m, n), one obtains, 
G fc:m ,„ » 8 fe — / #W ^e^ m+ ^« (l-(3/4)(0 2 + v9 2 -0v^)) 



27T i-00 27T J -00 

1 /-oo 1 /-oo 

— / #- 

27T J-00 27T j-00 

8 71 " -4(m 2 +w 2 +mre)/(9fc) 



i k — I rlrh— I dtp e i4,m+i,fin e -( 3fc / 4 )(^ 2 +v 2 -^) 



— 4(m^+n^+mn;/(9A;J (68") 



(2tt) 2 3^fc 

Substituting this expression into Eq. (65) and applying the Taylor expansion 
technique, identically as in Eq. (46), we obtain for the multiplicity distribu- 
tion of k gluons 

N (k) _ 8 fc 512 mn (m + n) ir f C F N c \ 

V " " {2-k? 81 V3^ eXP l C A k ^ ) m 

Eq. (69) is the final result of this sub-section. It has several interest- 
ing features. Firstly, unlike the case of many quarks, it has no dependence 
whatsoever on the cubic Casimir. Indeed, the result is entirely symmetric in 
m and n. Secondly, the weight in the exponential is just simply a Gaussian 
in the color charges. These are of order y/k ^> 1, and are therefore classical 
color charges. The prefactor, as in the SU(3) quark case, is proportional to 
mn(m + n). Therefore, as in the SU(3) case, one can write the sum over 
color charges in a box as an integral over an 8-dimensional classical mea- 
sure. We further note that for the gluons, the total number of degrees of 
freedom does turn out to be 8 fc even in our large k approximation. Finally, 
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the variance of the quadratic Casimir is proportional to an additional factor 
of Ca/Cf compared to the quark case. One thus obtains for the argument 
of the exponential (see Eq. (32)), 

C f Ca n2 c f r (Axi) 4 

c~ A T Q = c2 Ca ^T pp 

1 ^(Axjyy*, (70) 



2N C g 2 A 

where, in our simple model k is now the number of gluons in a tube of 
transverse size (Ax±) 2 , namely, k = (N 2 — 1) A (Axj_) 2 /irR 2 . One can thus 
deduce that the color charge squared per unit area of the gluon charges is 

V'A, g lue=9 2 N c —jp. (71) 

This result for the gluons is precisely the one expected in the MV model [28] . 
The rest of the derivation of the path integral for the MV-effective ac- 
tion goes through as discussed previously. We end this discussion with a 
caveat. It has long been understood that a renormalization group treatment 
where one includes bremsstrahlung gluons from quark sources, as additional 
sources for further bremsstrahlung of softer "wee" gluons, will modify the 
simple Gaussian of the MV model [20, 29, 21]. Naively, one can interpret 
this in terms of adding gluons to the color sources. The reason we have a 
Gaussian distribution of sources in the problem discussed here is two-fold; 
firstly, we ignore (quantum) correlations in the sources and secondly, the 
ground state in our problem is a singlet 1 representation. A closer analogy 
to the bremsstrahlung scenario would be to have quarks in the ground state 
and then add 8 representations. This would correspond to different initial 
conditions for our recursion relations and is outside the scope of this paper. 
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4.3 Many qq pairs 

In this case, we are multiplying 3x3 = 1 + 8 representations. Then, from 
the gluon case, it is straightforward to see that 

((1,0) x(0,l)) x (m,n) 
= (m, n) + (m + 1, n + 1) + (m + 2, n — 1) 
+ (m + 1, n - 2) + (m - 1, n + 2) 
+ (m, n) + (m, n) + (m — 2, n + 1) + (m — 1, n — 1) 

(72) 

The recursion relation for the multiplicity distribution N m>n is identical with 

the gluon case with the exception that one replaces 2N m>n — ► 3A^ mn in 

Eq. (56). The special cases are also very similar to the gluon case. The 

basic unit in the solution of Eq. (72) are the "n-nomial coefficients" in the 

expansion of the 9 k states as 9 k = (1 + 1 + 1 + 1 + 1 + 1 + 3) fe . 

The rest of the derivation is very similar to the previous cases and we 

will merely quote the result here. We obtain in the large k limit, 

at /i,n 9 k 27 \f?>mn (m + n) n m 2 +mn+» 2 

/vw ps ^ — e 2i (73) 

(2tt) 2 8 A; 4 l ' 

Again, one finds an expression symmetric in m and n, the ubiquitous factor 

mn(m + n) in the prefactor, and the Gaussian distribution of color charges. 

Since D™' n th {m 2 + n 2 +mn) /3, one can now easily deduce that fi\ = = 2 fi\. 



5 SU(N, 



We can carry out a similar analysis of the multiplicity of higher dimensional 
representations for the general SU(N C ) case. For problems in the standard 
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model, it is sufficient to have explicit solutions for SU(2) and SU(3) cases. 
However, there are several reasons to extend the analysis to the more general 
N c case. Firstly, we are presented with an interesting non-trivial random 
walk problem in N — 1 dimensions where the number of possible directions 
at each step depends on the type of representation being multiplied at that 
step. Secondly, we would like to confirm that the patterns we see in the 
N c = 2 and N c = 3 cases are not accidental but follow from a more general 
solution. For example, we interpreted the factors 2 and 3 that multiply the 
quadratic Casimirs in Eqs. (26) and (49) as the corresponding N c factors. 
We would like to verify in this section that they indeed are N c factors and 
not merely an accident of combinatorics. Another pattern we observed in 
previous sections is that the multiplicity is predominantly determined by 
the quadratic Casimir and that the most probable distribution has weights 
of order 0(l/y/k). We also would like to see that this is a general pattern. 
Finally, the large N c limit constitutes another classical limit of a spin system. 
Although this large N c limit is not what one needs in the small x MV model, 
working out an explicit expression for a general N c may shed some light on 
when exactly we can regard a system as a "classical" system. 

5.1 N c 

Let us consider adding a large number of quarks in the fundamental represen- 
tation of SU(N C ). For SU(N C ), a representation R can be uniquely labeled 
by an integer vector 

R=(r 1 ,r 2 ,---,r Nc -i). (74) 
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The corresponding Young tableau has N c rows with /j number of boxes in the 
i-th row. The number of boxes in the rows and r» are related by Vi = fi — fi+\ 
and /j's should satisfy the inequality 

fi>f2>---f Nc -i>f Nc (75) 

since t-j > 0. In this notation, a single quark is labeled by 

JV c =(l,0,..-,0) (76) 

We begin with a random representation given by 

R=(ri,r 2 ,---,rN c -i) (77) 

We would like to see what irreducible representations are generated when 
A/" c and R are multiplied together. 

Multiplying by N c corresponds to adding a box to the Young tableau for 
R. Adding a box to the i-ih row lengthens it by 1 or /$—>/$ + 1. Therefore 
the labels change as 

n = fi- fi+i ^r t + l (78) 

n-i = fi-i - /i ->• n - 1 (79) 

The result of the multiplication by N c is 

N c 

RxN c = J2 R 'm (80) 

m=l 

where 

R ' m = (ri,r 2 ,---,r Nc - 1 ) -e m _i + e m (81) 
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Here e m is a unit vector in the m-th direction (with e = ejv c = 0). As before, 
R' m with negative entries should be discarded. If one keeps multiplying iV c 's, 
the multiplicities of the representations in the k-th iteration and (k— l)-th 
iteration are related by 

N c 

G (k \R) = Y, G(k ~ 1) ( R + e rn^i-e m ) (82) 

m=l 

If k >> 1 and r^ ^> 1, we can make the continuum approximation by Taylor 
expanding the right hand side and discarding higher derivative terms: 

G {k) (R) « N c G (fc_1) (R) + b Nc G {k ^ ] (R) (83) 

where we defined 

N c -1 N c -2 

DN c =Y, d l-Y, d m d m+ i , (84) 

m=l m=l 

with d s = d/dr s . 

We also would like to make a similar continuum approximation for k as 
well so that Eq. (83) becomes a partial differential equation. Since N c can 
be large, care must be taken to ensure the validity of the Taylor expansion. 
To do so, we first let 

G ik \R) = G(k, R) = N k a{k, R) , (85) 

to take care of the potentially large iV c dependence. In the large k limit, the 
equation for a(k, R) is then 

d k a = j-rD Nc cr (86) 
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Letting the expressions we obtained for the SU{2) and S77(3) cases 
guide us, we try the following ansatz: 

a(k,R)=g(k)e W (-^p-\ , (87) 

where D 2 (R) is the quadratic Casimir of the representation R and r is a yet 
to be determined constant. Substituting Eq. (87) into Eq. (86) and solving 
for g(k) is tedious but straightforward. The solution is 



k (N c -i)/2 -r y 24k c k J ' 

where C is a normalization constant. Details of this derivation are in Ap- 
pendix C. 

Note the appearance of N c D 2 jk in the exponential. Both the SU(2) and 
SU(3) cases we worked out match up with this expression and the expression 
for the color charge squared per unit area obtained in section 3.2 



2nR 2 ' 



^ = ^1' f8()) 



is valid for any N c . The additional term in the exponential only affects the 
overall normalization and is negligible when N c <C k as in the £77(2) and 
577(3) cases. 

The expression in Eq. (88) is a particular solution of the recursion rela- 
tion in the large k limit. However, it is not guaranteed that Eq. (88) satisfies 
the necessary boundary condition. In fact, the solution (88) corresponds to 
the Gk-.m,n functions in the previous section. The multiplicity N(k, R) should 
therefore be a particular linear combination of G(k, R) as in Eq. (41). 

Unfortunately, we have been unable to find the analog of Eq. (41) for 
general N c although it is conceivable that a general solution may be obtained 
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along lines presented in Ref . [38] and Ref . [39] . However, if the pattern seen in 
SU(2) and SU(3) were to persist, one might guess that in the large k limit, 
the multiplicity should be given by 

N(k, R) = Afd(R) exp (-N c ^^\ (90) 

where M is a normalization constant and d(R) is the dimension of the R 
representation in SU(N C ). We can also expect the color phase space measure 
will have the structure 

c^c-iq = dndr2 . . . rfr ^_ 1 d(Rf dVt Nc (J v c _ 1} (91) 

where dQ]y c (jv c -i) ar e the R independent canonically conjugate Darboux vari- 
ables. We have been unable thus far to verify these conjectures explicitly. 

5.2 N c x N c 

We can easily generalize the above analysis to the case of a large number 
of quark and anti-quark pairs. For SU(N C ), the fundamental and its dual 
representations are labeled by 

JV C = (1,0,...,0) (92) 

JV C = (0,0,--.,1) (93) 

As previously, we start with a representation 

i*=(ri,r 2 ,...,rx_i) (94) 

Multiplying by iV c results in 

( r i> r 2> • • • » r N c -i) = ( r i> r 2, ■ ■ ■ , r Nc -i) ~ e m _i + e m (95) 
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for < m < N c . Further multiplying by N c results in 

l r l ? r 2> ' ' ' ) r Af c -l/ = VI' r 2' ' ' ' ) r N c ~l) "I" e n-l — e n (9o) 

for < n < N c . 

Hence the multiplicities in the fc-th iteration and the (k— l)-th iteration 
are related by 

G«(i*) = g G'( fe - 1 Hi2 + e m _ 1 -e m -e n _ 1 + e n ) (97) 

m,n=l 

In the continuum limit, where all the t-j are typically large (or k is large), 

/JV C -1 W c -2 \ 

G^(i2) = iV c 2 G( fe - 1 )( J R)+2iV c £ ^,- £ 3L«Ui CC" 1 '^) (98) 

\ m=l m=l / 

We again let 

G (*-i)(jj) = G(A;, i?) = JV?* 0(jfe) exp f-M^P) (99) 

\ rk J 

The solution for the resulting differential equation is 

Again, this is a particular solution of the recursion relation. The comments 
at the end of the last section (after Eq. (88)) are equally applicable here. 

5.3 The adjoint representation 

The adjoint case isn't much different from the quark- ant iquark case. For 
SU(N C ), 

(iV 2 -l) = (l,0,---,0,l). (101) 
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We know that 

N c x N c = 1 © (iV c 2 - 1) (102) 

So all we have to do is to change Eq. (97) to 

G^ k \R)= £ G'( fe - 1 )( J R + e m _ 1 -e m -e n _ 1 + e n )-G'( fe - 1 )(i2) (103) 

m,n=\ 

and the solution is 

Again, this is a particular solution of the recursion relation. The comments 
after Eq. (88) also apply here. 

6 Discussion and Summary 

We have discussed in this paper the distribution of color charge representa- 
tions generated by k partons in an SU(N C ) gauge theory. Since the partons 
are random, the problem is a random walk problem in the space spanned by 
the SU(N C ) Casimirs. We explicitly considered the SU(2) and SU(3) cases 
before considering the general N c limit. For all the cases considered, we find 
that the most likely representation is one of order 0{yk). The distribution 
of representations about this representation is given by an exponential in the 
quadratic Casimir D2, with a weight proportional to k. In the case of SU(3) 
quarks, the contribution due to the cubic Casimir is suppressed relative to 
the leading Gaussian term by 0{l/yk). Remarkably, for gluons and quark- 
anti-quark pairs, the result is given entirely in terms of the quadratic Casimir. 
Our results for SU(3) can be generalized to the SU(N C ) case. Although the 
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random walk is now in the space of N c — 1 Casimirs, the distribution of 
representations is still given by the quadratic Casimir. 

Since the most likely representation is of order 0{yk), one can argue 
that, in the yk — > oo limit, the representation is a classical representation. 
The arguments here are no different from those for the treatment of classical 
spins. A formal representation of SU(N C ) representations, suitable for path 
integral formulations, can be made in terms of coherent states [32]. As dis- 
cussed in appendix A, it has been shown in Ref. [31] that, for large k, one 
recovers the classical limit. 

Our discussion here was formulated in the framework of the McLerran- 
Venugopalan model for small x physics. The Gaussian distribution of color 
charge sources and their treatment as classical color charges was assumed 
in this model. We have provided here a more rigorous basis for these as- 
sumptions. The small x behavior of QCD is more complex than the one 
outlined in the MV model and more sophisticated renormalization group 
(RG) treatments have been developed [14, 21]. A key feature of the MV 
model, namely, the lack of correlations among the sources, breaks down in 
these treatments. Therefore, unsurprisingly in our random walk picture, the 
distributions are no longer Gaussian. The classical assumption persists how- 
ever since the most probable representations are likely still classical ones. 
This assumption deserves further study. Remarkably, in the limit of very 
small x, a mean field treatment of the RG equations, recovers a Gaussian 
distribution of sources-albeit a non-local one [12, 13]. Interactions among 
sources in this limit produces screened charges of size 1/Q S , where Q s is 



37 



the saturation scale 13 . In our language, this might mean that the effective 
charges are random and can therefore be represented by Gaussian sources. 
It is interesting to speculate whether the combinatorial techniques developed 
here for higher dimensional representations can be used in numerical studies 
of high energy Onium-Onium scattering in QCD [15]. 

Other clear applications of the techniques developed here are in trans- 
port problems in QCD. It has been argued previously that the transport of 
color in high temperature QCD could be represented in terms of the classical 
dynamics of color charges [27]. It was however not clear when this classical 
description was the appropriate one and attempts to derive these from first 
principles have had to resort to ad hoc approximations [36]. We do know 
though that any systematic treatment of transport problems involves coarse 
graining of color charges over large distance scales. These can therefore be 
treated using the recursion techniques developed here. These developments 
will be discussed elsewhere [37]. 
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A Large k color charge representations as clas- 
sical color charge representations 

An interesting limit of quantum theories is the limit of N — > oo, where N 
can denote either the underlying invariance group of the theory or the size of 
higher dimensional representations when the invariance group of the theory 
is held fixed. In both cases, the large N limit corresponds to a classical limit 
and can be formally shown to correspond to the limit where h — ► 0. As 
discussed in a nice review by Yaffe [2], the quantum dynamics of a system 
will reduce to classical dynamics as h — »■ 0, if and only if the system can be 
prepared in a state whose uncertainty in its conjugate momenta and positions 
vanishes in this limit. 

Coherent states [30] have this property and it is therefore useful to write 
states of the quantum theory in this representation. These coherent states 
are orbits of the Heisenberg-Weyl group, and can therefore be generalized to 
any Lie group. For an SU(N C ) gauge theory, one can define coherent states 
carrying SU(N C ) color charges [32]. 

In Ref. [31], coherent states for SU(N C ) groups were studied with the 
particular purpose of studying the classical limit. The authors identify co- 
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herent states of a particular representation in terms of the signature of rep- 
resentation (which is none other than size of the totally symmetric tensor-in 
our case k [33]) and the weight vectors n of the representation. They define 
the variance of the square of the length of the isospin vector (the quadratic 
Casimir D 2 = J2Ql), 

AD 2 = (tt| £Qal*> - ]T(M>|Q a |M>> 2 = (*p2 - D' 2 \m) , (105) 



In the coherent state basis for fixed SU(N C 

AD 2 N c 



(106) 



D 2 N c + k 

This result is valid for k ^> iV c . Thus the large k limit, where the variance 
vanishes can be identified with the h — > limit in quantum mechanics, and 
one can formally relate h = 1/k. We refer interested readers to Ref. [31] for 
a more extensive discussion of this correspondence. 

B The volume of the canonical phase space 
in SU(3) 

In Refs.[26, 27], the phase space variables for 577(3) are parameterized as 

Ql = COS0i7T + 7T_, Q 2 = sin0 1 7T + 7T_ 
Q?> = Hi 

Q 4 = C++-K+A + C + -n-B, Q 5 = S ++ n + A + S + -n_B 
Q 6 = C-+H-A-C—1I+B, Q 7 = S- + 7r-A-S—7r+B 
Qs = tt 2 (107) 
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with 

and 
A = 



COS ( -(±01 + V3(j)2 ± 03) ) , <S±± = Sin ( -(±01 + V3(j)2 ± 03 



(10* 



2vr 3 \ 



/ J I — J 2 7T 2 \ 

I 3 + - + V3J 


I 3 + *> + js) 


/2 Ji + J 2 
V 3 


-7T 3 - 


7T 2 \ 



5 



1 
2^\ 



J 2 -Jl 7T 2 WJi + 2J 2 7T 2 W2Ji + J 2 7T 2 \ 

(109) 



where (Ji, J 2 ) equals the (m,n) label of a SU(3) representation. The pairs 
(iTi, 0i) form the canonical pairs. The total phase space volume given for a 
fixed (m, n) is given by the integration over these canonical variables: 



Q(m,n) 



> i d7ii d<p2 dit2 d<p3 dir 3 



110) 



Evaluating the integral over the conjugate momenta is not trivial. It can 
be inferred from Fig. 1 of Johnson's paper [26] but the explicit derivation, as 
shown below, is quite involved (if straightforward). 

The values of the canonical momenta 7Tj are restricted to make A and 
B real. To solve for the ranges of 7r's, it is convenient to define 



K x 



2m+n 7v- _ 2n+m 



X = 7T 3 + f 3 , y = 7T 3 -f 3 



'MD 
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which also gives 



2 
dxdy = d(Tc 3 + rc 2 /V3)d('K 3 — 7r 2 /v3) = —j=d^^d^ 2 (H2) 

The factors in A and B are 

A ± = K x -K 2 + x (113) 

A 2 = K 2 + x (114) 

A 3 = K x -x (115) 

B 1 = K 2 -K x +y (116) 

B 2 = K 2 -y (117) 

£? 3 = ^i+y (H8) 

Potentially, there are 16 sign combinations for the factors A t and Bi that 
make the products A\A 2 A 3 and B\B 2 B 3 positive. For instance, to make the 
argument of A positive, one may require all the factors in A to be positive 
or require that A\ and A 2 to be negative and A 3 to be positive. 

Fortunately, due to the conditions x + y > 0, K x > 0, K 2 > and 

2K X - K 2 = m > (119) 

2K 2 - K x = n > (120) 

one can easily show that the only consistent sign combination is for all At 
and Bi to be positive. 

The condition that all Ai are positive yields the range 

K 2 - K x < x < K x (121) 
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The condition that all Bi are positive yields the range 

K x -K 2 <y<K 2 (122) 

then 

dx dy dni = dx dy (x + y) 

K-i-Kx JK1-K2 J-(x+y)/2 JK2-K1 JK 1 -K 2 

mnjm + n) , . 



Equivalently, 



a/3 
d7rid7r2d7T3 = ——mn(m + n) (124) 



The integrals over the angles can be evaluated in the following way. We 
first note that 



,2 V 



J d 8 Q5{Q 2 - 1) = fdQQ 7 f dQ 7 5{Q 2 

(125) 



7^ 
6 ' 



where we have used the fact that the area of a unit sphere in 8-D is 27r 4 /3! 14 . 
On the other hand, using Eq.(51) this can be also written as 

£ = Jd*Q6(Q*-l) 

= / d(p\d(p2d<pzdn\dn2dii3dmdnmn('m + n)— — 5 1 

J 48 y 3 I 

(126) 

We already know the result of the 7Tj integration from Eq. (124). Therefore, 

7r 4 1 [,,,,,, 22/ 2 ( m 2 + mn + n 2 \ 

— = — / dcpidfadcpzdmdnm n (m + n) 1 (127) 

6 64 j V 3 ) 



14 For a nice discussion of the volume and area of SU(N) groups, we refer the readers 

to [34]. 
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The delta function gives 



mi = (128) 



and the condition m > becomes n < y/3. Hence 



,, 2 2 ( m 2 + inn + n 2 \ f^s 3(m 2 n 2 (m + n) 2 ) 



dmdn m n (m + n) 6 1 = / dn 



o 2m + n 

^ , VSn 2 (-3 + n 2 f 
an 



m=m\ 



which results in 

2(2tt) 



>HJup2U,ipz 



o V4 — n 2 

2?r . , 

Tr (129) 



(130) 



V3 ' 

The total phase space volume defined in Eq. (110) is therefore determined 
from Eq. (124) and Eq. (130) to be 

(27r) 3 
fi(m,n) = mn(m + n). (131) 

C Young tableaux and the quadratic Casimir 

in SU(N C ) 

In this section, we establish the relationship between a Young tableau and 
the quadratic Casimir used in the text to obtain SU(N C ) results. 

For SU(N C ), a Young tableau has N c rows. Suppose a Young tableau has 
fj boxes in the j'-th row. Then this tableau corresponds to a representation 
labeled by 

fl=(ri,r 2 ,--.,r JVc _i) (132) 
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where r,i = f\ — /j+i. Inverting this relationship yields 



with 



fi = T,r jt (133) 



1 / JV <=- 1 \ 

r Nc = — {K- £ jrJ , (134) 

where 

^ = E/* ( 135 ) 

is the total number of boxes in the tableau. 

According to Okubo [35], the quadratic Casimir for SU(N C ) is given by 

D 2 (R) = ^ (136) 

where 

N c 

C 2 (R) = K(l + N c -K/N c ) + ^f j (f j -2j) (137) 

i=i 

Since D 2 is a Casimir that depends only on R, it must be independent of K 

even if the expression (137) explicitly contains K. The fact that C 2 is in fact 

independent of K can be easily checked by taking a derivative with respect 

to K and using that fact that only r^ depends on K. 

In the main text, the equation we need to solve has the form 

d k a(k,R) = -^D Nc a(k,R) (138) 

where 

N c -1 N c -2 



m=l m=l 



DN c ^J2 d l-J2 d mdrn + l (139) 

7? 
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Our ansatz is 



a(k,R)=g(k)e X p(-^^\ (140) 

The left hand side is then 

d k a = d k g e~ c ^ k + ^g e~ c ^ rk (141) 

where we used the fact that C 2 is independent of K or equivalently, k. The 
right hand side is 



-C 2 /rk = n jj^ p -C 2 /rk 
1 



D Nc ge-w = gD Nc e~ 

fN c -\ N c -2 \ 

E (<^ 2 ) 2 - E (d m C 2 )(d m+1 C 2 ) ge~ c ^ rk 



\m=l m=l / 

1 /JVc-1 N c -2 \ 

- -T E 9 mC 2 - E ^9 m+1 C 2 ^e"^ (142) 

rfc \m=l m=l / 

Since the expression Eq. (137) depends on /j, we need to know 

d 
dmfi = 7; — (rj + r l+1 H h r^J 

= (1 - m/iV c ) - 0(Z - m) (143) 

where 9(1 — m) = 1 if / > m and 9(1 — m) — if I < m. Therefore, 

d m C2 = 2 E fidmfi — 2 E ^m/j 

iV~c iV c 

= 2g/,((l- m/7V c ) - 0(Z - m)) - 2 E * ((1 - W^c) - 0(1 - m)) 
1=1 1=1 

m 

= -2(m/N c )K + m(N c -m)+2j2fi (144) 

1=1 

where we used the fact that X^=i fi — K. We note that this formula works 
even if m = N c or m = 0. In that case 8qC2 = 0, 0n c C2 = 0. 
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Second derivatives are 

m 

d 2 m C 2 = d m 2j2fi = 2d m (K- E fi. 



N c 



1=1 

-2 E (-- 

l=m+l ^ ^ c 

m(N c — m) 



l=m+l 



(145) 



and 



N c 



d m +id m C 2 



-2d m+ \ 2_^ fi 

l=m+l 

m(N c — m — 1) 



(146) 



Using these expressions, it is tedious but straightforward to work out 

N c -1 N c -2 

E &A - E d m+1 d m C 2 = (N c - 1) (147) 



ra=l 



ra=l 



and 



N c -1 



N c -2 



Y(d m C 2 ) 2 - Y, ( d m+iC2)(d m C 2 ) = 2C 2 + ^^- ^ (148) 

m=l m=l 



The equation (138) then becomes 

ag 



n ^ ° 2 

d k g + -r^g v 

rk z N r 



r 2 k 2 l 6 rk 



2a 
N~ r . 



d k \ng 



N C 2 (N C 2 - 1) 1 



We first let 

then solve 

which yields 

where S is an arbitrary normalization constant. 



2Aak 2 



#-': 



' ! " ' ' 2ak 



(149) 



(150) 



(151) 



(152) 
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